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Conventional Cooper pairing arises from attractive interaction of electrons in the metallic bands. 
Recent experiment on Co-doped LiFeAs shows superconductivity in the insulating valence band, 
which is evolved from a metallic hole band upon doping. Here we examine this phenomenon 
by studying superconductivity in a three-orbital Hamiltonian relevant to the doped LiFeAs. We 
show explicitly that Cooper pairing of the insulating hole band requires a finite pairing interaction 
strength. For strong coupling, the superconductivity in the hole band is robust against the sink of 
the hole band below the Fermi level. Our theory predicts a substantial upward shift of the chemical 
potential in the superconducting transition for Co-doped LiFeAs. 


In conventional BCS theory, superconductivity arises 
as a Fermi surface instability, and superconducting (SC) 
gap decreases exponentially with the density of states at 
the Fermi energjffl. So fully occupied or empty bands 
are usually ignored in study of superconductivity. It is 
interesting that the opening of a SC gap on a fully occu¬ 
pied band is observed in a recent ARPES experiment in 
a Fe-based superconductoil^. Fe-based superconductors 
have generated great interest because they have high¬ 
est transition temperature next to the cuprates at am¬ 
bient pressur^^ltiSl, Typically, a Fe-based superconductor 
has a complex Fermi surface consisting of both electron 
pockets and hole pocketJ^^^^. In comparison with other 
iron pnictides, LiFeAs has a much shallower hole pocket 
around the center of the Brillouin zon^i^lb^. Upon Co¬ 
doping, the entire hole band sinks below the Fermi en¬ 
ergy, the hole pocket disappears and the hole band be¬ 
comes insulating as revealed in the ARPES experiment. 
Below the SC transition temperature Tc, the energy gap 
between the top of the hole band and the Fermi energy 
is found to be larger than the gap in the insulating state, 
which suggests Cooper pairiim and superconductivity in 
the fully occupied hole band .12 

It is interesting to note that a similar phenomenon was 
discussed theoretically by Nozieres and PistolesP^ over 
a decade ago. They considered superconductivity in a 
semiconductor with a small band gap, and found that the 
semiconductor may undergo a transition into SC state at 
low temperature if the binding energy of the Cooper pair 
is larger than the energy cost to produce free charge car¬ 
riers across the semiconducting gap. LiFei_a;Coa; is a 
multi-band metal with electron band metallic and hole 
band fully occupied. The physics behind the supercon¬ 
ductivity of the fully occupied valence or hole band is 
similar to the superconductivity in a semiconductor. In 
the multi-band Fe-based compound, electrons in the fully 
occupied hole band may be excited to the states in metal¬ 
lic electron bands, which introduces an additional pairing 
channel. This allows the realization of the band insulator 
to superconductor transition in multi-band LiFei_a;Coa;. 
In this paper, we consider a phenomenological mean field 
Hamiltonian based on a three-orbital model, relevant to 


LiFei-xCoa;, and study the Cooper pairing and supercon¬ 
ductivity of the model. We show explicitly that Cooper 
pairing in the insulating hole band requires a finite pair¬ 
ing interaction strength. For strong pairing interaction, 
the superconductivity in the hole band is robust against 
the sink of the hole band below the Fermi level. Our the¬ 
ory predicts a substantial upward shift of the chemical 
potential in the SC transition for Co-doped LiFeAs. 

We consider the Hamiltonian 

TL ='Ho + Hpair, ( 1 ) 

where Hq and Hpair describes the kinetic part and the 
attractive interaction respectively. In the LiFeAs mate¬ 
rial, the states around the Fermi energy mainly consist 
of dxz, dyz, and dxy orbitals of Fe-3d orbitals. So we 
consider a three-orbital model introduced by Brydon et 

al.l^ 

’?^o(k) = T (k)C]{^g.,^CkCTn — M cj^.^^Ckcrm (2) 

kcrmn kcrm 

where cj^.j^^(cko-m) creates (annihilates) an electron with 
orbital m, crystal wave-vector k and spin a, m,n = 1, 2, 3 
corresponds to dxz, dyz and dxy orbital respectively. The 
hopping integrals are given by 

= 2ti cos ky + 2^2 cos kx + cos kx cos ky 
+ 2tii {cos2kx — cos2fcj^) 

= 2ti cos kx + 2^2 cos ky + cos kx cos ky 
— 2tii (cos 2kx — cos 2ky) 

= A^xy + 2^5 (cos kx + cos ky) + dte cos kx cos ky 
+ 2tg (cos 2kx + cos 2ky) 

+ 4tio (cos 2kx cos ky + cos kx cos 2ky) 

= 4^4 sin kx sin ky 

= 2 it 7 sin kx + ditg sin kx cos ky 

= 2itg sin ky + ditg sin ky cos kx (3) 

where ti are the hopping constants and chosen to be the 
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same as in the paper by Brydon et aP^. The band disper¬ 
sions and the Fermi surfaces for LiFeAs are depicted in 
Fig-B There are three bands, the a band with smaller 
hole Fermi pocket, the /3 band with larger hole Fermi 
pocket, and the 7 band with electron Fermi pockets. The 
chemical potential in the normal state is determined by 
the relative position of the top of the a band. We choose 
/i = 0.358met^ to model LiFeAs in the absence of Co¬ 
doping. In this case, the chemical potential is 2meV be¬ 
low the top of a band, consistent with the AREPS data. 
The partial substitution of Co by Li introduces electrons 
into the FeAs plane, which c an be treated effectively as 
lifting the chemical potentiaP^^IIlIl jn our calculations 
below, we choose fi = 0.364meE and — O.SGSmeE to 
model the Co-doped cases in experiments with the chem¬ 
ical potential AmeV (1 per cent of Co-doping) and 8meV 
(3 per cent of Co-doping) above the top of the a band, 
respectively. 
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FIG. 2. SC and insulating phase diagram of the valence 
a band at T = 0, in parameter space Vi and V 2 for a fixed 
AE = AmeV, (a); and in parameter space of Vi and AE for 
a fixed V 1 /V 2 = 1.2, (b). Vi and V2 are intra-band pairing 
coupling and inter-band pair hopping, respectively, and AE 
is the Q-band gap, namely the energy difference between the 
top of the a-band and the chemical potential p. AE = 4meV 
corresponds to Co-doping x = 0.01. In these calculations, fj, 
is fixed. 
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FIG. 1. (a): Band structure obtained from the tight binding 
model for LiFeAs. Red, green, and blue curves are for a, /3, 
and 7 band respectively. Black dashed line in the inset indi¬ 
cates the chemical potential in the absence of Co-doping, (b): 
Fermi surfaces in the model with two hole pockets around F 
point and two electron pockets around X and Y point. In 
the calculations, the parameters in Eqn. (2) and (3) are cho¬ 
sen as, in units of eV, = 0.02, t2 = 0.12, fa = 0.02, t4 = 
-0.046, ts = 0.2, fe = 0.3, ty = -0.15, fs = = 

-0.06, Co = -0.03, til = 0.014, = 1, /r = 0.358. 


the superconductivity. We assume that the attractive 
interaction is from the next nearest neighbor intersite 
pairin^^. Then the gap function in k-space has the form 
cos kx cos ky , and the pairing Hamiltonian reads 

■Hpair = - ^ cos kx cos ky -I- h.c.^ , (4) 

k,r=Q!,/3 

where A'^ is the amplitude of the spin singlet SC gap in 
T = a,j band, which are defined as 

A“ = ^ ^ cos k^ cos ky (Vi (c£^c“k;) + E2 (cilt^-kf)) 

k 

A^ = ^ ^ cos cos ky (Vi + V ^2 (c£^c“k;)) • 

k 

( 5 ) 


The self-consistent equations then read 


The band structure calculated in the tight binding 
model is in good agreement with the ARPES results. 
However, we note that there is a degeneracy of the a 
and 13 bands at the F point in our model, which is absent 
in the ARPES data. Since the Co-doping mainly affects 
the a band but not much about the /3 band as observed in 
the ARPES, we shall drop out the /3 band for simplicity 
in the calculations below to study the superconductivity 
of the insulating a band. Namely, we will consider only 
the a-hole band and the 7 - electron band. The role of the 
/3 band will be briefly discussed at the end of the paper 
and it is similar to the 7 band. 

We now consider superconductivity. Eor this purpose, 
we introduce a phenomenological pairing term Hpair- Be¬ 
cause the a band and 7 band do not have any overlap 
near the Fermi energy, we may safely neglect the inter¬ 
band pairings, and include the intra-band pairing cou¬ 
pling and the inter-band pair hopping process to study 
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i ^ J pQ:(7) 


X tanh 


El^^''l2kBT 


( 6 ) 


where 


0(7) 

M 


[A“(t') cos kx cos ky 


and are the quasiparticle energy and single particle 
energy of 0 ( 7 ) band, respectively. 

We first discuss the zero temperature (T = 0) case. 
The superconductivity in the a band is controlled by Vi, 
V 2 j and the a band gap AE. For simplicity, we choose a 
gauge where A“ is real and non-negative. By solving the 
self-consistent equation for various Vi and V 2 at a fixed 
AE = 4meV, we find a phase diagram plotted in Fig. 
da). The SC gap in the a band, A“ increases with the 
increase of the coupling constant, and Vi and V 2 play a 
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FIG. 4. Chemical potential shift 5^ in the SC state as 
function of the insulating a band gap Ai? (panel a); and 
the SC gaps on the a band with 5^ included (red dashed 
curve) and neglected (blue solid curve), and the SC gap on 
7 -band (green dashed line), (panel b). The parameters are 
l/i = 76meV and = 1.2 


to an upward shift of the chemical potential as shown in 
fig- n The chemical potential shift may be calculated 
through the equation for number of electrons per site 


FIG. 3. Schematic illustration of the electronic bands, panels 
(a) and (b), and the mean occupation number of the electrons 
of a state, panel (c) and (d), in the normal and SC states. In 
the SC transition, some of electron pairs in the fully occupied 
valence band transfer to the conduction band, which leads to 
Cooper pairing in the insulating hole band. Chemical poten¬ 
tial in the normal state is shifted upward to fi + in the 
SC state due to the transferred electrons into the conduction 
band. 


similar role in enhancing the superconductivity. The de¬ 
pendence of the A“ on AE and Vi is plotted in Fig. [^b) 
for a fixed value of V 1 /V 2 = 1.2 . From Fig. 3, it is clear 
that the insulating a-band becomes SC only at Vi > I 4 
with Vc the critical pairing coupling. This conclusion is 
similar to the study on semiconductor pairing^i^. A“ re¬ 
duces as Vi reduces to and vanishes when Vi < 14. In 
conventional BCS theory, the superconductivity instabil¬ 
ity is induced by infinitesimal coupling constant. For an 
insulating valence band, one must remove some electrons 
from the fully occupied valence band to the conduction 
band, which will cost an energy at least equal to the insu¬ 
lating band gap AE. Therefore, the a band may become 
SC only when the energy gain by forming a Cooper pair 
is larger than the energy cost, which leads to a finite 
critical Fe¬ 
in the above calculations, the chemical potential /r is 
fixed. In the conventional BCS theory for metallic nor¬ 
mal state, the chemical potential shift is tiny if we fix 
the total electron density. The situation is different in 
the LiFei_ 2 ;COa;As case. As we analyzed above, some 
electrons must be removed from hole band to induce su¬ 
perconductivity in the a band. If we consider the total 
electron density is fixed as in the usual case, the only 
place those electron can go is the 7 band. That increases 
the number of electrons in the electron band and leads 
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+ 



(7) 


We have solved eqn. ® and eqn. Q for the couplings 
Vi = 76meV and V 1 JV 2 = 1-2 at various a band gap 
AE. The qualitative behavior on the superconductivity 
is similar to that we discussed at the fixed chemical po¬ 
tential, and the SC gap reduces with increasing a band 
gap as shown in fig. |^b). For the present model, as AE 
increases from 4 meV to 8 meV, A“ reduces from 6 meV 
to 1.8 meV. The calculated chemical potential shift is 
depicted in fig. |^a), which is rather large especially in 
the case with larger A“. For example, 6^ ~ l.lmeF with 
AE = 2meV, while the SC gap on a band is just around 
3.9meV. Therefore, the chemical potential shift can not 
be simply neglected in this case. Similar with the SC gap, 
the chemical potential shift also reduces with increasing 
the a band gap AE. The reason is that the total number 
of electrons excited from a band to 7 band reduces with 
suppression of the SC gap, which leads to the reduction 
of Sfi. 

The chemical potential shift may be important for 
the estimation of the SC gap from the experimental data. 
If one neglects J/r, the SC gap on the a band can be 
calculated from the gap at F point in the normal state 
(AE) and gap in the SC state (Aafter sc), 


Ao = V^A^ftecsc-Ai?^ (8) 

Including Sfj, , the SC gap is given by 

A“ = v/A2^tecsc-(Ai5 + ^M)^ (9) 

From the above analyses, the SC gap estimated without 
including the chemical potential shift is larger than the 
correct one. In fig. |^b), we depict the value of Ag (blue 
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solid line) and the correct SC gap A“ (red dashed line) 
for comparison. 



FIG. 5. Temperature dependence of SC gap on a and 7 
band from model calculations. In upper panel, Vi = 76meV 
, V 1 IV 2 = 1.2. (a): AS = 2meV\ and (b): AS = 8 meF. In 
lower panel, Vi = IQOmeV, V 2 = ImeV . (c): AS = 2meV ; 
(d): AS = 8meV. 

We now study superconductivity at finite tempera¬ 
tures. The SC gap is reduced with the increase of tem¬ 
perature as shown in fig. where the results for various 
Vi, V 2 , and AS are depicted. As shown in the figure, be¬ 
cause of the presence of interband pair hopping, the SC 
gap of a band and 7 band vanishes at the same temper¬ 
ature Tc though they are different at zero temperature. 
If one sets V 2 = 0, the a band and 7 band are decoupled 


in the SC transition, and they may have different Tc- 

In the above calculations, we have neglected the /3 hole 
band. The /3 band will play a similar role as the electron 
band to the a band superconductivity. The /3 band will 
enhance the superconductivity on the a band through 
the pair hopping between the two bands. The effect to 
the chemical potential shift is mixed. The enhancement 
of superconductivity will increase the number of electrons 
in the a band to be excited, hence to increase the chem¬ 
ical potential shift upward. On the other hand, the f3 
band provides additional states for the a electrons to 
transfer to, which will increase the density of states at 
Fermi energy and reduce the chemical potential shift. So 
the quantitative result depends on the parameters. But 
the qualitative result will not change. 

In summary, we have presented a theory to explain the 
observed superconductivity in the a-hole band which is 
completely below the Fermi energy upon Co-doping in 
the LiFei_a;Coa;As. We consider a three-orbital Hamilto¬ 
nian with an intra-band spin-singlet pairing and an inter¬ 
band pair hopping terms. The a band is shown to be¬ 
come superconducting if and only if the pairing strength 
exceeds a critical value. Therefore the observation of 
superconductivity in an insulating valence band is an ex¬ 
perimental proof of the strong coupling in Fe-based su¬ 
perconductors. For strong pairing interaction, the super¬ 
conductivity in the hole band is robust against the band 
gap of the hole band. We show that the chemical poten¬ 
tial shift accompanied with the superconducting transi¬ 
tion is rather large and should be taken into account in 
extracting the superconducting gap from the ARPES ex¬ 
periment. 
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